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I. INTRODUCTION 

Defining a Yang-Mills theory at the conformal 3 + 1 
spacetime boundary of an anti-de Sitter (AdS 5 x S 5 ) 
space [1] demands massless mesons. But it was shown 
that nonperturbative calculations can be performed in 
gauge theories with a mass gap dual to supergravity in 
warped spacetimes [2]. One can truncate the AdS space 
(hard-wall) or introduce a potential in the AdS space 
(soft-wall) to break the conformal invariance to build 
confinement at large distances, while retaining confor¬ 
mal behavior at short distances. The mass spectrum of 
hadrons was computed by Teramond and Brodsky [3] by 
truncating the AdS space (in its holographic coordinate 
z) with a hard-wall set by the QCD scale (A qcd)- 

The correspondence between light-front holographic 
QCD and AdS§ space was further utilized when hadronic 
wave functions were obtained by mapping an impact 
variable Q which represents the measure of transverse 
separation of the constituents within the hadrons, to the 
holographic variable z [4]. Four-dimensional Schrodinger 
equations were obtained for the bound states of massless 
quarks and form factors were evaluated for spacelike Q 2 . 

The hard-wall model is the simplest but not the best 
way to incorporate confinement. The conformal invari¬ 
ance can be broken by introducing a dilaton background 
in the theory. It was found that by choosing a specific 
profile for a nonconstant dilaton, the usual Regge de¬ 
pendence can be obtained [5]. The dilaton field <p(z) 
produces an effective confining potential in the AdS ac¬ 
tion. Choosing <p(z) = n 2 z 2 , usual Regge behavior was 
observed and meson spectrum and properties have been 
studied within holographic QCD frameworks by several 
groups [6-9]. AdS/QCD wave functions and distribution 
amplitudes have been used to explain rho meson electro¬ 
production [10] and radiative B decays [11]. 

The AdS/QCD model consists of massless quarks due 
to lack of a corresponding field in AdS action which could 
generate quark masses. A prescription was suggested by 
Brodsky and Teramond [12] to include quark masses, 


which was further developed into models to obtain me¬ 
son wave function with massive quarks [13, 14]. In this 
work, the model developed in Ref. [14], which incorpo¬ 
rates massive quarks, is used to look at wave functions 
and distribution amplitudes for mesons. Then, form fac¬ 
tors and transition form factors are computed with mas¬ 
sive quarks and compared with experimental data. The 
model is also studied for massless quarks and compared 
with the massive one. 

The paper is organized as follows. Original AdS/QCD 
model for massless quarks is discussed briefly in Sec. 
II. Then, massive quarks are introduced and developed 
upon in Sec. II A. Wave functions and distribution am¬ 
plitudes are studied in Sec. Ill for 7 r, p , K, and J/if 
mesons. Then we calculate form factors for the four 
mesons and charge radii for pion and kaon and compare 
the form factors and charge radii for pion and kaon with 
the experimental data in Sec. IV. Further, transition 
form factors are computed for photon-to-meson decay of 
7r, 77 , and rf and plotted against experimental data in 
Sec. V. Finally, y 2 per degree of freedom is calculated 
to quantitatively compare the predictions of the massive 
and the massless quark models in Sec. VI. The work is 
concluded in Sec. VII. 


II. MESON IN ADS/QCD 

Working in the Fock basis, the meson wave equation 
is obtained within semiclassical approximation with all 
the interaction terms embedded in an effective potential 
U(Q. The light-front wave equation for the meson is 
then obtained as [ 6 ] 

^+ u (o)m= M2 m, a) 

where ( is a light-front transverse variable specifying 
the separation between the quark and the antiquark. 
Whereas, the wave equation in AdS space for meson of 
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spin- J is given by 


, 3—2 j 


a <p{z) 


oV{z 


, 3-2 J^z 


(■ mR) z 


<M-) = M 2 <f>j(z). 

( 2 ) 

Factoring out the dilaton factor from the AdS field as 
/ d \ J—3/2 

*j{z) = (jJ e~^ 2 Mz), (3) 


and identifying z —> £, we get the following form for 

u(0-. 

u (0 = ^"(0 + |^(C ) 2 + (4) 

Using the dilaton profile of the soft-wall model <p(z) = 
k 2 z 2 , the Schrodinger equation for the meson is obtained 
as 

(-$- 1 1c^ + kV + 2k2(j - 1) ) So 

= M 2 0(C), (5) 


with the mass spectrum 

M 2 = 4k 2 ^n + J + ^ , (6) 


where n and L are the radial and orbital quantum num¬ 
bers. 

The effective light-front wave function (LFWF) for a 
two-parton ground state in impact space comes out as 

<?W(C) ~ -\A(1 - x)e~^ K ' c \ (7) 

In Eq. (7), using £ 2 = x(l — x)b 2 ± to replace £, where 
b_L is the transverse impact variable that is conjugate to 
the light-front relative transverse momentum coordinate 
kj_, the meson LFWF can be rewritten as (along with 
some constants chosen for this model [14]) 

b±) = ^ y /x(l-x)e(-^ 1 -^), (8) 

V7T 


with a normalization parameter A, which will be fixed 
later. Performing a Fourier transform, we get 


</W(u kj_) = 


4nA 


Ky/x( 1 — x) 




(9) 


In a completely different approach, namely, a domain 
model of QCD [15], meson spectra and wave functions 
have similarities with the soft-wall AdS/QCD. It will be 
interesting to map these two formalisms to have a better 
understanding of the origin of the confining potential or 
the particular choice of the dilaton profile. 


A. Massive quark model 


The quark masses (mi and m 2 ) are introduced by 
extending the kinetic energy of massless quarks with 


Kn = 


x(l—x) 


to 


„ „ „ O Q 7/4-1 

Kq —> K = Kq + /i 12 , H 12 = 


m 2 
1 — x 


This is equivalent to the change 

d 2 d 2 2 
~d( 2 “dC 2 +Ml2 


Finally the LFWF comes out to be 

/ f v ^ 4lrA 

V , qi q2(X,k±) = 


K\J x(l — X) 


x exp — 


2k 2 x{1 — x) 


2k 2 J 


.( 10 ) 


This modified wave function, which incorporates quark 
masses within the soft-wall framework, will be used to 
investigate the properties of ir, p, K , and J/if) mesons. 


B. Getting the parameters 


The dilaton parameter, k, is fixed by the Regge tra¬ 
jectory. Now that quark masses have been introduced, 
the mass spectrum in Eq. (6) will be modified too. The 
modified spectrum should look like [16] 


M 2 = 4k 2 




xf 2 (x, mi, m 2 ) + other corrections, 


( 11 ) 


where /(a;,mi, m 2 ) is the wave function for the longitu¬ 
dinal mode, which is factored out from the meson LFWF 
as 


= -^^=e lMv f(x,m 1 ,m 2 ). (12) 

V 2ttC 

For massless quarks, f(x) has been found to be f(x) = 
\/x(l — x ), which follows from an analysis of the meson 
form factor [4] and is extended for massive quarks as 


f(x,m 1 ,m 2 ) = Nf(x)e ^ 


where N is normalized as 1 = dxf 2 (x,m\,m 2 ). This 

and other corrections to the mass spectrum are discussed 
in detail in Ref. [17]. Other corrections come from one- 
gluon exchange and hyperfine-splitting contributions to 
























3 


the effective meson potential £/(£). However, eventually 
after calculating all the corrections, it is seen that all the 
additional terms do not require a change in k to fit the 
experimental spectrum for the case of light mesons [17]. 

For 7 r, p, and K , we take n = 540 MeV [18] and for 
J/ip, n = 894 MeV [14], which are obtained by fitting the 
meson mass to its respective Regge trajectory. Since 7r, 
p, and K follow a Regge trajectory different from that 
followed by J/ip, they have different values of k. There 
are some other values of k for the pion in the literature, 
e.g., n = 375 MeV [6] and 361 MeV [13], obtained by 
fitting the pion form factor, and k = 432 MeV [19], ob¬ 
tained by fitting the pion transition form factor. Here, 
we consider n = 540 MeV, obtained through the Regge 
trajectory and it is the same for 7r, p, and I\. 

Masses and decay constants are taken from PDG and 
other sources. To fix the normalization parameter A, we 
use the expression for the decay constant [20] 

f M =2V&J dx j ^^V> 9l q 2 (a:,k ± ). (13) 

All the values of the parameters used in this work are 
listed in Table I. The pion has its parameters listed 
for zero, current, and constituent quark masses, respec¬ 
tively, and the kaon has its parameters listed for current 
and constituent quark masses, respectively, while p and 
J/ip have their parameters listed only for constituent 
quark masses. 


Meson 

mi, m 2 (MeV) 

k (MeV) 

A 

/m (MeV) 

7T (u,d) 

0 , 0 

540 

0.79 

130.7 


5, 5 

540 

0.79 

130.7 


330, 330 

540 

2.24 

130.7 

P (u,d) 

330, 330 

540 

2.65 

154.7 

K (u,s) 

5, 95 

540 

0.99 

156.1 


330, 500 

540 

4.69 

156.1 

J/-0 (cc) 

1500 

894 

684.4 

277.6 


TABLE I. Parameters used in this work. 


III. WAVE FUNCTION AND DISTRIBUTION 
AMPLITUDE 

For 7r, p, K , and J/t/h the wave functions are plotted 
in Figs. 1-6 using Eq. (10). These graphs show the 
wave function narrowing down along x as the mass of 
the meson increases. This is because of the decreasing 
momentum spread among the quarks for heavier mesons. 
For current quark masses, the wave functions peak near 
the end points (x = 0,1), while for constituent quark 



FIG. 1. The pion wave function ^(r, ki), for current quark 
masses mi = m 2 = 5 MeV. 
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FIG. 2. The pion wave function ip(x, kj_), for constituent 
quark masses mi = m 2 = 330 MeV. 


masses, which are heavier compared to current quark 
masses, the wave functions peak near x = 1/2. Also, 
Figs. 4-5 for the kaon wave function show an asymmetry 
in the momentum distribution because of unequal quark 
masses. 

The meson distribution amplitudes are important for 
theoretical description of electroproductions of mesons. 
The distribution amplitude is calculated as [21] 

ri 2 d 2 k± 

4‘{x,q)= 0u k ±), (14) 

where q 2 is the transferred momentum squared. Al- 
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FIG. 3. The p wave function ^>(x, kj_), for constituent quark FIG. 5. The kaon wave function kj_), for constituent 

masses mi = m 2 = 330 MeV. quark masses m\ — 330, m 2 = 500 MeV. 
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FIG. 4. The kaon wave function ^(x, kj_), for current quark 
masses mi = 5, m 2 — 95 MeV. 


FIG. 6. The J/^> wave function kj_), for constituent 
quark masses mi = m 2 = 1500 MeV. 


though the meson 1^) can be expanded into Fock states 
I'*/’} = ai\qq) + a 2 \qqg) + a 3 \qqgg) H-, for asymptoti¬ 

cally large q 2 , the first term dominates and we can use 
the wave function in Eq. (10), which is for the qq state, 
and integrate over kj_ to get <j>(x) = <j>(x, Q oo) as 

■Kx) = ^ Vx(l - x) exp (- x) ) . (15) 

The distribution amplitudes for 7r, p , K , and J/'ip are 
plotted in Figs. 7-12 in which they are compared with 
the pQCD prediction [22] q oo) = y/3/2fMx(l — 
x). Here, again, as the meson gets more massive the 


distribution narrows down toward the center. 


IV. FORM FACTOR AND CHARGE RADIUS 

In light-front coordinates, the electromagnetic form 
factor can be computed from the matrix elements of 
the plus-component of the current J + at light-front time 
x + = 0 as 


(P , |J + (0)|P) = (P + P') + F(q 2 ), 


(16) 
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FIG. 7. The pion distribution amplitude for current 

quark masses, as calculated from the model (solid line) and 
from the pQCD prediction (dashed line). 



FIG. 8. The pion distribution amplitude 4>{x), for constituent 
quark masses, as calculated from the model (solid line) and 
from the pQCD prediction (dashed line). 



FIG. 9. The p distribution amplitude for constituent 

quark masses, as calculated from the model (solid line) and 
from the pQCD prediction (dashed line). 



FIG. 10. The kaon distribution amplitude for current 

quark masses, as calculated from the model (solid line) and 
from the pQCD prediction (dashed line). 


where q 2 < 0 is the transferred spacelike momentum 
squared. J + {x ) can be computed as 

J + {x) = ^2e q tjj(x)^ + ilj{x). (17) 

<2 

Using the momentum expansion of 4’+(. x ) = A +ip(x ) = 
7 ° 7 + 2 /;(a;) in terms of creation and annihilation operators 


[23] 

i>+(x~,x±) a = J ^r[ 6 A(9)wa(g,A)e“* 9 - x 

+d\(q)^v a (q, X)e lq ' x ], (18) 

J + {x) can be expressed in the particle number repre¬ 
sentation. To compute the matrix element in Eq. (16), 
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Then, after using a normalization condition and inte¬ 
grating over intermediate variables in the q + = 0 frame, 
we obtain the Drell-Yan-West expression [24, 25] for the 
form factor 


-Fm(<7 2 ) =J2 [ ^ [d 2k -L 

n j 

X 'tpn/M {p^i o , A^) • (20) 

A Fourier transform gives the form factor in terms of q 
as 


n— 1 

Fm(<1 2 ) = ZhU 

n j= 1 


dxjd 2 b±j exp 


'j 


^ 'jn/M (Xj: b -Lj ) i 


( 21 ) 


FIG. 11. The kaon distribution amplitude rf)(x), for con¬ 
stituent quark masses, as calculated from the model (solid 
line) and from the pQCD prediction (dashed line). 


where j runs over all n constituent quarks. 

For a meson, n = 2 and there are two terms which 
contribute to Eq. (21). So, we get 



-Fm(<7 2 ) [ dxd2h - L e’ qj -' a:b - L |i/’„/M(a;.b±)| 2 

n,j 

= [dxd 2 b ± e^ xcose \ip n/M (x,b ± )\ 2 . 

n,j 

( 22 ) 


Now, the integral over b_L gives a Bessel function of the 
first kind and the sum over two quarks runs for charges 
e i | and for momentum fractions x —1 x, 1 — x 

giving two integrals. Exchanging 1-iin the second 
integral we get 


F M (q 2 ) = 2tt [ 
Jo 


dx 


x(l — x) 
x [ (d( J 0 ( C<7 



\^m{x, C )| 2 , 


FIG. 12. The 3/ip distribution amplitude </>(*), for con¬ 
stituent quark masses, as calculated from the model (solid 
line) and from the pQCD prediction (dashed line). 


the initial and final meson states |V , m(-F >+ > P_l)) are ex¬ 
panded in terms of their Fock components 



x |n : XiP + , XiPx + kj_i, 


Vhi. (xi- kj_f, Aj) 

A»). (19) 


(23) 

where ( 2 = x(l — x)b^. 

Now, performing a Fourier transform on our modified 
LFWF in Eq. (10) and using ( 2 = x(l — x)bj_, we get 
LFWF in impact space including light-quark masses as 


ip M ( x , C) = 47tAk^/x(1 



2 


(24) 


Using Eq. (24) in Eq. (23), we get the meson form factor 
as 
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F M (q 2 ) = (327t 3 AV) 



,- k 2 r 


(25) 


A. 7r form factor 

For the pion, the form factor F^Q 2 ) ( Q 2 = —q 2 > 0) 
is plotted in Fig. 13 for constituent quark masses. To 
compare with the massless quark model, the pion form 
factor with zero-quark masses is also plotted. Here, 
Fn(Q 2 ) is normalized to F w (0) = 1. The pion form fac¬ 
tor has been measured by CEA/Cornell [26], CERN [27], 
Baldini et al. [28], JLAB [29-31], and CLEO [32, 33]. 
The data from Baldini et al. [28] have many points which 
overlap with the ones from CEA/Cornell [26] and CERN 
[27] and, hence, have been plotted only once. We com¬ 
pare our model with these experimental data sets up to 
Q 2 = 18 GeV 2 in Fig. 13. Clearly, the pion form factor 
with our model fits better with the data than the one 
with zero-quark mass limit. For large values of Q 2 , the 
asymptotic pQCD prediction for the meson form factor 
is [34] 

Fm(Q 2 ) = 8lTa ^ + . (26) 

For a s = 0.3 and f M + = f n + = 130.7 MeV, the pion 
pQCD data are shown in Fig. 13. The result from 
our model with massive quarks agrees well with the 
asymptotic pQCD prediction. This agreement is better 
than that with the data for large Q 2 values from CLEO 
[32, 33], 

B. p form factor 

For p, the form factor F p (Q 2 ) is plotted in Fig. 14, 
along with the plot for zero quark masses. Here, F p {Q 2 ) 
is normalized to F p ( 0) = 1. 



FIG. 13. The pion form factor FV(<3 2 ) for nonzero (solid 
line) and zero (dashed line) quark masses with k = 540 MeV 
[18]. The plot magnified for smaller values of Q 2 is shown 
in the inset. The data are taken from [26-33]. The data 
points plotted from Baldini et al. [28] (shown as black aster¬ 
isks) are only those which do not overlap with the ones from 
CEA/Cornell [26] and CERN [27]. 


D. J/if> form factor 

For J/ij), the form factor Fj/^(Q 2 ) is plotted in Fig. 16, 
for both nonzero and zero quark masses. Here, Fj^(Q 2 ) 
is normalized to Fj/^,{ 0) = 1. 

E. Charge radius 


C. K form factor 


The charge radius of the meson can be calculated as 


For the kaon, the form factor Fk(Q 2 ) is plotted in 
Fig. 15 for both nonzero and zero quark masses. Here, 
F k (Q 2 ) is normalized to Fk( 0) = 1. The kaon form 
factor has been measured by CERN [35], FERMILAB 
[36], and CLEO [32, 33]. We compare our model with 
these experimental data sets up to Q 2 = 18 GeV 2 in 
Fig. 15. Similar to the pion, the values from the asymp¬ 
totic pQCD prediction [Eq. (26)] are also plotted, which 
fit closely with our model. 


„ 2 \ 


6 

m 


dF(Q 2 ) 

dQ 2 


Q 2 —0 


(27) 


For the pion, with n = 540 MeV, the charge radius comes 
out to be a/ (r 2 ) = 0.529 fm for massive quarks and 1.679 
fm for massless quarks, while the experimental value is 
y/ir 2 ) = 0.672 ± 0.008 fm [37], 

For the kaon, with k = 894 MeV, it comes as a/ (r 2 K ) = 
0.598 fm for massive quarks and 1.679 fm for massless 



























FIG. 14. The p form factor F P (Q 2 ) for nonzero (solid line) 
and zero (dashed line) quark masses. 



FIG. 16. The J /'ip form factor Fj/^(Q 2 ) for nonzero (solid 
line) and zero (dashed line) quark masses. 



limit. In this range, our model does not agree quite well 
with the experimental data, which in turn causes this 
mismatch. 

V. TRANSITION FORM FACTOR 

The amplitude for two-photon production of the me¬ 
son ( 7*7 —I M) contains an unknown function of pho¬ 
ton virtuality ( Q 2 ) which is known as the meson transi¬ 
tion form factor. The photon-to-meson form factors are 
measured in several experiments [38-42]. Thus, there 
are many theoretical interests to predict these transition 
form factors. Here we present the predictions from the 
AdS/QCD model of massive quarks. 


FIG. 15. The kaon form factor Fk(Q 2 ) for nonzero (solid 
line) and zero (dashed line) quark masses with k = 894 MeV 
[14], The plot magnified for smaller values of Q 2 is shown in 
the inset. The data are taken from [32, 33, 35, 36]. 


quarks, while the experimental value is \J(r^) = 0.560± 
0.031 fm [37]. 

The deviation from the experimental values, signif¬ 
icantly evident in the pion’s case, can be attributed to 
the fact that the value of a charge radius mainly depends 
on the values of the form factor only at small Q 2 (near 
zero), given that we compute the slope in the Q 2 = 0 


A. 7r transition form factor 


The photon-to-pion transition form factor F n ^(Q 2 ) 
can be calculated in QCD as [21] ( Q 2 = — q 2 > 0) 


Q 2 F ni (Q 2 ) 


4 

VS Jo x 


1 + 0 



(28) 


where x is the longitudinal momentum fraction of the 
quark hit by the virtual photon in the hard scattering 
process and x = 1 — x is for the other spectator quark. 
The distribution amplitude <f>(x, Q ) is given by Eq. (14). 
In pQCD, the pion distribution amplitude for Q 2 —> 00 
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is given by cj)(x,Q 2 —> oo) = ^/3/2f 7I x( 1 — x ). Thus, This simplifies to the following form 
the pQCD prediction for the pion transition form fac¬ 
tor is given by lim q^^ 00 Q 2 F^{Q 2 ) = where f n 

is the pion decay constant. The Q 2 dependence of the 
form factor helps to constrain the models for distribution 
amplitudes. In our model, the distribution amplitude is 
given by Eq. (15) and the expression for the pion tran¬ 
sition form factor can be written as 



The pion transition form factor is plotted in Figs. 17-18 
for F V1 {Q 2 ) and Q 2 F^(Q 2 ), both for nonzero and zero 
quark masses. The pion transition form factor has been 
measured by BABAR [38], CLEO [39], CELLO [40], and 
Belle [41]. We compare our model with these experimen¬ 
tal data sets up to Q 2 — 40 GeV 2 in Figs. 17-18. The 
plot for our model with nonzero quark masses fits bet¬ 
ter with the data. The data for large Q 2 values from the 
BABAR Collaboration do not agree so well with the the¬ 
oretical prediction from our model, similar to what has 
been seen in several other studies. However, the data 
from the Belle Collaboration for the similar range of Q 2 
values fit better with the model. Note that the form fac¬ 
tor with massless quarks fails to fit with the experimental 
data. 

B. r/ and r/ transition form factors 

The form factors for rj and rj mesons, for the photon- 
to-meson decay, can be calculated using the prescription 
used in Ref. [19]. The r\ and rf mesons result from the 
mixing of neutral states rjg, and 771 given by 

uu + dd — 2 ss uu + dd+ ss 

718 = w ; Vl = 71 ' 

Their transition form factors are given by the same ex¬ 
pression as for the pion, except there is an overall con¬ 
stant factor cp = ^ and ^(2 for the r)s and 771 mesons, 
respectively. By multiplying cp with Eq. (30), the rjs 



FIG. 17. The pion transition form factor F JT1 (Q 2 ) for nonzero 
(solid line) and zero (dashed line) quark masses. The data 
are taken from [38-41]. 


and 771 transition form factors are obtained which give 
the transition form factors for the physical states 77 and 
77 ' through the following mixing of eigenstates 

( F V1 \ l cos 9 — sin0 
F r] t 1 J l sin0 cos 9 


F„ 


( 31 ) 
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FIG. 18. The pion transition form factor Q 2 F JT1 (Q 2 ) for FIG. 19. The r/ transition form factor F V ^(Q 2 ) for nonzero 

nonzero (solid line) and zero (dashed line) quark masses. The (solid line) and zero (dashed line) quark masses. The data 

data are taken from [38-41]. are taken from [39, 42]. 


The mixing angle is 9 = —11.4° [37]. The plots for 
FM-yiQ 2 ) and Q 2 F Ml (Q 2 ) for rj and rf are shown in 
Figs. 19-20 and Figs. 21-22, respectively. These transi¬ 
tion form factors have been measured by BABAR [42] 
and CLEO [39]. We compare our model with these ex¬ 
perimental data sets up to Q 2 = 40 GeV 2 in Figs. 19-22. 
The data from the BABAR Collaboration for rj and rf 
have better agreement with the model’s prediction than 
in the case of 7r. For the transition form factors, the 
quark masses seem to play a major role and, particu¬ 
larly at large Q 2 , significant improvements are observed 
over the model with massless quarks. 


VI. CALCULATION OF y 2 PER DEGREE OF 
FREEDOM 

Before looking at the accuracy of the predictions of 
this model, it needs to be emphasized that the purpose 
of this work is to demonstrate the improvement in agree¬ 
ment with the experimental data of the AdS/QCD model 
for meson, when the wave functions are modified to in¬ 
clude the quark masses. The model still represents a 
valence quark picture and, hence, is not expected to re¬ 
produce the experimental data with high accuracy. From 
Figs. 13-22, it is evident that the inclusion of quark 
masses considerably improves the model’s predictions. 
The only tunable parameter in the AdS/QCD model is 
k, which is already fixed by the Regge trajectory, and, 
hence, we have no free parameter to fit. For the sake of a 



FIG. 20. The r/ transition form factor Q 2 F v ~ t {Q 2 ) for nonzero 
(solid line) and zero (dashed line) quark masses. The data 
are taken from [39, 42]. 


quantitative comparison, in Table II, we have listed the 
values of y 2 per degree of freedom for different meson 
properties that we have computed in this paper. Since 
we have compared two models, the uncertainty in the 
value of k, which is the same for the massive and the 
massless quark models, is not considered in the y 2 com¬ 
putation. 
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FIG. 21. The rj transition form factor F v i^(Q 2 ) for nonzero 
(solid line) and zero (dashed line) quark masses. The data 
are taken from [39, 42]. 



FIG. 22. The 7 / transition form factor Q 2 F rl i 1 (Q 2 ) for 
nonzero (solid line) and zero (dashed line) quark masses. The 
data are taken from [39, 42]. 


From Table II it can be inferred that, except for the 
pion and the kaon form factors, the massive quark model 
gives much better values of x 2 P er degree of freedom 
for the meson properties. The massive quark model re¬ 
produces the overall behavior of the pion form factor 
much better than the massless one, as can be seen from 
Fig. 13. However, as shown in the inset of Fig. 13, for 


very small values of Q 2 , where we have a large cluster of 
data points, the massive quark model does not fit well 
with the data and, hence, x 2 /d.o.f. becomes larger for 
it than that for the massless case. For the kaon form 
factor, in the inset of Fig. 15, we have magnified the 
small Q 2 region, which clearly shows the improved re¬ 
sult of the massive quark model over the massless one. 
Here, again, x 2 /d.o.f. comes out to be large, primarily 
due to three data points with very tiny errors at large 
Q 2 values, where the massive quark model deviates con¬ 
siderably. When we remove these three points from the 
analysis (remember that we are not fitting any parameter 
here), the x 2 /d.o.f. value for the massive quark model 
improves drastically over that for the massless one, as 
shown in the Table II. 


Meson 

property 

X 2 /d.o.f. 
(massive quarks) 

X 2 /d.o.f. 

(massless quarks) 

F n (Q 2 ) 

73.3715 

21.5423 

Fk(Q 2 ) 

174.4685 

23.333 

F k {Q 2 ) a 

0.819 

4.4018 

f.AQ 2 ) 

6.6594 

104.4687 

FvAQ 2 ) 

2.1186 

62.8514 

Frt'AQ 2 ) 

4.0627 

327.4731 


a Here, we exclude the three data points at large Q 2 values 
shown in Fig. 15 

TABLE II. Value of y 2 per degree of freedom for different 
properties of mesons presented in Figs. 13-22. 


VII. CONCLUSIONS 

Light-front holographic QCD maps an effective grav¬ 
ity theory defined in a five-dimensional warped anti- 
de Sitter spacetime to a semiclassical approximation to 
strongly coupled QCD quantized on the light-front. It 
provides a way to model the wave functions of hadrons, 
analyze their spectrum, and their properties and com¬ 
pare with the experimental data. The soft-wall holo¬ 
graphic model has been found particularly useful in com¬ 
puting hadron form factors and transition form factors, 
which are found to be consistent with the currently avail¬ 
able experimental data. 

In this paper, we have presented the results for the 
wave functions, distribution amplitudes, and form fac¬ 
tors for 7r, p , I \, and J/ip mesons using the wave func¬ 
tions predicted by the soft-wall holographic model in 
AdS/QCD. Here, we have used wave functions, modi¬ 
fied to include quark masses, which follow from the pre¬ 
scription suggested by Brodsky and de Teramond [12]. 
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The results for pion and kaon form factors are compared 
with the available experimental data from different ex¬ 
periments. In most cases, AdS/QCD results with zero 
quark mass fail to agree with the experimental data. The 
modified AdS/QCD wave functions for massive quarks 
considerably improve the agreement with the data. The 
distribution amplitudes and form factors are also com¬ 
pared with pQCD predictions and found to be consis¬ 
tent. In fact, the pQCD predictions for the pion and 
kaon form factors, given by Fm(Q 2 ) = 87 ra s f^+fQ 2 
at large Q 2 values, are in excellent agreement with our 
AdS/QCD model with massive quarks. 

We have also presented the photon-to-meson transi¬ 
tion form factors for 7 r, 77 , and 7 /, which are reasonably 
consistent with the available experimental data. The 
BABAR data [38] at large Q 2 values for the pion tran¬ 
sition form factor are not compatible with our model’s 
results, as has been the case with many other theoreti¬ 
cal models, while the data from Belle Collaboration [41] 
have better agreement. The BABAR data [42] for rj and 
7 / are more compatible with the model’s results than 
they are for n. Again, the results for massless quarks 
are inconsistent with the data. These conclusions are 


also supported by the values of y 2 per degree of free¬ 
dom, computed for the massive and the massless quark 
models. 

Overall, the results for the form factors and the tran¬ 
sition form factors predicted from our model with mas¬ 
sive quarks are seen to be more consistent with the ex¬ 
perimental data than the results from the model with 
massless quarks. Inclusion of quark masses can be seen 
as a good approach to extend the already available 
AdS/QCD models in order to explain hadronic proper¬ 
ties, e.g., generalized parton distributions, nucleon form 
factors, etc. 
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